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A system of coupled scalar fields is introduced which possesses a spectrum of 

(N ; 

O , massive single-soliton solutions. Some of these solutions are unstable and decay 

' into lower mass stable solitons. Some properties of the solutions are obtained using 

C ■ 

^ ' general principles including conservations of energy and topological charges. Rest 

>—i ; 

\^ , energies are calculated via a variational scheme, and the dynamics of the coupled 

fields are obtained by solving the field equations numerically. 

(N . 
> ■ 
(N ■ 

(N : 

(N ■ I. INTRODUCTION 

o ■ 
(N : 

^ \ Relativistic solitons, including those of the conventional sine-Gordon equation, exhibit 

' remarkable similarities with classical particles. They exert short range forces on each other 
Oh! ^ 

O . and make collisions, without losing their identities They are localized objects and do 

^ \ not disperse while propagating in the medium. Because of their wave nature, they do tunnel 
• 1—1 , 

^ ■ a barrier in certain cases, although this tunnelling is different from the well-known quantum 

H ; 

_cd \ version . Topological solitons are stable, due to the boundary conditions at spatial infinity. 

Their existence, therefore, is essentially dependent on the presence of degenerate vacua |@. 

Topology provides an elegant way of classifying solitons in various sectors according to 
the mappings between the degenerate vacua of the field and the points at spatial infinity. 
For the sine-Gordon system in 1+1 dimensions, these mapping are between cj) = 2mT, n G Z 
and X = ±oo, which correspond to kinks and anti-kinks of the SG system. More complicated 
mappings occur in solitons in higher dimensions. For example, in cosmic strings, the vacuum 
is and topological sectors corresponds to distinct mappings between this and a large 
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circle around the string. This leads to the homotopy group 7r^(5'^) = Z. 

Coupled systems of scalar fields have been investigated by many authors ||, ^ . Bazeia 
et al. considered a system of two coupled real scalar fields with a particular self-interaction 
potential such that the static solutions are derivable from first order coupled differential 
equations. Riazi et al. ||^ employed the same method to investigate the stability of the 
single soliton solutions of a particular system of this type. Inspired by the well-known 
properties of sine-Gordon equation, we introduce a coupled system of two real scalar fields 
which shows a considerably reacher structure and dynamics. The present system is not of 
the form investigated in [Q and and static solutions are not derivable from first order 
differential equations. 

In our proposed system, a spectrum of solitons with different rest energies exist which 
are stable, unstable or meta-stable, depending on their energies and boundary conditions. 
Some of the more massive sohtons decay spontaneously into stable ones which subsequently 
leave the interaction area. Note that the term "soliton" is used throughout this paper for 
localized solutions. The problem of integrability of the model is not addressed here. Such 
non-dispersive solutions are called "lumps" by Coleman ||^ to avoid confusion with true 
solitons of integrable models. However, it has now become popular to use the term soliton 
in its general sense. 

Pogosian |Q] investigated kink solutions in bi- dimensional SU{N) x Z2 models. He found 
that heavier kinks tend to break up into lighter ones. Comparing our results with those 
reported by Pogosian, it is interesting to note that fairly similar phenomena are observed in 
quite different systems. In an earlier paper, Pogosian and Vachaspati reported (A^ + l)/2 
distinct classes of kink solutions in an SU (N) x Z2 field theory. 

The organization of the paper is as follows: In section II, we introduce the lagrangian 
density, dynamical equations, and conserved currents of the proposed model. In section 
PI , some exact solutions, together with the corresponding charges and energies are derived. 
The necessary nomenclature, and general behavior of the solutions due to the boundary 
conditions are also introduced in this section. Numerical solutions corresponding to differ- 
ent boundary conditions are presented, and properties of these solutions like their charges, 
masses and their stability status are addressed in this section. In order to investigate the 
stability of the numerical solutions, their evolution are worked out numerically. The dy- 
namical evolution of the unstable solutions are investigated further in section |V[ Our final 
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conclusion and a summary of the results is given in the last section. 



II. DYNAMICAL EQUATIONS AND CONSERVED CURRENTS 

Our choice of the Lagrangian density reads 

£ = ^^'^019^01 - ai(02' + ei)(l - cos0i) + ^9^029^02 - a2(0i' + e2)(l - cos02) (1) 

in which a\^i and ei^2 are positive constants, and 0i^2 are two real scalar fields. The back- 
ground space-time is assumed to have the metric g^y = diag(l, —1) in 1+1 dimensions and 
c = 1 has been used throughout this paper. Recall that the sine-Gordon Lagrangian density 

£ = -dydf,(f) - a{l - cos 0) (2) 

leads to single soliton solutions 

= 2n7r±4tan-^e^^ (3) 

all having rest energies 8^/a. Our proposed Lagrangian density comes from the idea 
of mixing two scalar fields in such a way that the discrete vacua at 0i = 2mr, 02 = 2m7r 
(m, n G Z) survive, while the rest energies of solitons for each field (say 0i) is affected by 
the value of the other (02) and vice versa. It will be seen later, that this idea leads to the 
appearance of non-degenerate solitons with distinct topological charges. 

The corresponding equations of motion are easily obtained by applying the variational 
principle 6{ J Cd?x) = to the lagrangian density (|]): 

□01 = ai(02^ + ci) sin 01 + 2a20i(l - cos02), (4) 

and 

□02 = 02(01^ + £2) sin 02 + 2q;i02(1 - COS0i). (5) 
Since the lagrangian density (0) is Lorentz invariant, the corresponding energy-momentum 



tensor IC 



Tf"" = dyid''(j)i + 9'^029>2 - g^'^C, (6) 
satisfies the conservation law 

d^T"" = 0. (7) 



FIG. 1: The self-interaction potential is shown as a height diagram over the (01,02) plane. 



The Hamiltonian density is obtained from Eq. (|^) according to 

n = T^= i0l' + ^02' + ^0'i' + ^0'2' + V^(01, 02), (8) 

where dot and prime denote time and space derivatives, respectively. It is seen from (|l]) 
that the potential is 

r(0i, 02) = ai(02' + ei)(l - cos0i) + a2(0i' + e2)(l - cos 02). (9) 

This potential is shown in Fig. (|1|), as a height diagram over the (0i,02) plane. In this 
figure and any numerical calculation throughout this paper, we fix the parameters as ai = 
0.3, a2 = 1, ei = 0.5, and €2 = 0.07 for the sake of being definite. The classical vacua 
consist of points in the (0i, 02) plane at which the condition V{(f)i, 02) = holds. Since the 
potential consists of two non-negative terms, it vanishes if and only if the two terms vanish 
simultaneously. This leads to the vacuum set of points in the (0i, 02) space: 

V = {(01, 02)101 = 2m7r and 02 = 2n7r; m,nGZ}. (10) 

It can be shown that the following topological currents can be defined, which are conserved 
independently, and lead to quantized charges: 

J':av = 5m,e^'^5.02/27r. (11) 
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In these equations, m is an integer, and the integers p and q are defined according to 

p = integer part of (0i/27r) + 1, 
q = integer part of (02/27r) + 1. 

The subscripts V and H denote "Vertical" and "Horizontal" which will be explained later. 
The currents J^jj y are conserved, independent of each other: 

d.J^v = 0. (12) 
The corresponding topological charges are given by 

/oo 
J^nd^ = 5mp [0i(+oo) - 0i(-oo)] /27r, 
-oo 

/oo 
J^ydx = [02(+oo) - 02(-oo)] /27r. (13) 
-oo 

These charges quantify the mappings between the vacua 0i_2(±C)o) G V and the points at 
spatial infinity. 

III. SINGLE SOLITON SOLUTIONS 

Static solutions which correspond to transitions between adjacent vacua are symbolically 
shown in Fig. (^. Accordingly, we call the static solutions H (horizontal) and V (vertical) 
types. We have called them 'horizontal' and 'vertical' simply because of their orientation in 
(01)02) plane. Initial guesses for these solutions are obtained using the known properties 
of the conventional sine-Gordon equation 0. Note that the exact solutions are not strictly 
horizontal or vertical lines in the (0i,02) plane. Rather, they are bent curves joining two 
neighboring vacuum points due to the coupling between the two scalar fields. Any finite 
energy solution should start and end at one of the vacuum points belonging to V (Eq. |10[). 
Equations (H) and (|^) possess the following exact single soliton solutions: 

Hp_^_i : 01 = 4tan~"'^ exp [±a(a; — Xq)] + 2pn and 02 = 0, (14) 

with 

QmH ^mpi QmV 0) 
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FIG. 2: Nomenclature of horizontal (H) and vertical (V) solutions according to boundary condi- 
tions. The three lightest solutions: Hf, Vq, and Hi are stable, while and V2 are unstable. The 
corresponding rest energies and charges are shown in Table | for our choice of parameters. 



and 



with 



^2 = 4 tan ^ exp [±b{x — Xq)] + 2q7i and 0i = 0, 



(15) 



QmH — 0, 



QmV — 



mqi 



where p and q are integers. Note that Eqs. (|T^ and (|T5|) satisfy Eqs. (||) and (||) in the 



static 



dt 



0) case. In these equations, a = ^a\e\, b = ^0262, and Xq is the kink position. 
Apart from these exact solutions, we will introduce other static solutions which are obtained 



numerically later. Despite the similarity of the special solutions (|Tj) and ([T5|) with those 
of sine-Gordon equation, there are profound differences between the general static solutions 
of the present system and SG, including non-degenerate soliton masses and instability of 
some of the static solutions. The rest energies of the H and y-type solitons are obtained by 
integrating the corresponding Hamiltonian densities (with 0i = 02 = 0) over the x-space. 
The rest energies are approximately given by: 



H -type: Mh ^ 8^/ai{A7v^n^ + ei), 



(16) 



and 



V - type : My ~ 8^/a2i^7i^n^ + €2), 



(17) 



where n is an integer. We have written programs in the 'MATLAB' environment, which 
calculate static solutions and the dynamical evolutions from pre-specified initial conditions. 
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FIG. 3: Examples of static solutions. 
The static solutions are obtained by minimizing the energy functional 

E = I Hdx (18) 



where 7i = Tq is the energy (Hamiltonian) density, using a variational procedure. The 
program is asked to terminate as soon as the energy converges up to a 10~^ accuracy. 
Inspired by the exact solutions and ([T5|) , the initial guesses were taken to be of the 
form of the sine-Gordon solitons. These initial guesses were deformed toward the lowest 
energy solutions by the variational calculation. This deformation is caused by the slope of 
the potential function as depicted in Fig. (|^). Sample solutions are shown in Fig. @. 
Static solutions Vq, Hf, Hi and as projected on the (0i,02) plane are shown in Fig. (H). 
A sample classification is shown in Table | for a particular choice of the parameters. Anti- 
solitons exist with the same masses (rest energies) as solitons but with opposite charges. 
The anti-soliton solutions are obtained by simply exchanging the boundary conditions at 
X = +00 into those at x = —oo. For example, 

^0 . J M-'^) = 0, 0i(+oo) = 27r ^^^^ 
^ b2{-oo) = 0, 02(+oo) = 0, 




7C 27t 



FIG. 4: The static solutions Vq*^, H^, and as projected on the {(f)i,<p2) plane. 
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TABLE I: A sample classification of the lowest energy solitons for ai = 0.3, a2 = 1, ei 
62 = 0.07. 
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while 

^0 . I M-^) = 27r, 0i(+oo) = ^^^^ 
^ I 02(-oo) = 0, M+oo) = 0. 

It is clear from the definition of topological charges (|Tl|) that solitons and anti-solitons 
have opposite charges. The approximate formulae (0) and (0) may be compared with the 
numerical results depicted in Table |. It can be seen that the masses of the H-tjpe solitons 
are better approximated by Eq. (|l^) than those given by Eq. (|T^) . This difference can be 
attributed to the fact that the horizontal solitons deviate less from the sine-Gordon solitons 
(see Fig.(l)). 



IV. SOLITON DECAY 



The time-dependent solutions of the dynamical equations (^) and (^) were obtained nu- 
merically, by transforming these PDEs into finite difference equations, and calculating the 
fields t) and </)2(x,t) in successive time steps. Most of the static solutions obtained did 
not undergo appreciable variations when inserted into our dynamical program as initial con- 
ditions. This is a numerical indication for the stability of the corresponding static solutions. 
However, not all the static solutions obtained by the variational method explained in the 
last section were found to be stable. For the given choice of parameters, V^, for example, is 
unstable and decays spontaneously via (see Fig. |(a) and (b)) 

V^' Hi + + Hi (21) 

Here, H^ is the anti-soliton of if J*. Figure (^a)) shows the 0i and 02 fields after several 
time steps. The decay product solitons can be identified at positions where the fields jump 
from one vacuum to an adjacent vacuum. Figure (^ helps recognizing these decay products. 
Figure (H) shows the stable solutions H\, Vq, and H^ and the unstable solution on the 
(01,02) plane. Numerical calculations show that the decay of starts with the trajectory 
shown for this solution and evolves gradually to the H^ [H^ in the reverse direction), V^^ 
and H\ trajectories in this figure. Figure (^(b)) shows the corresponding energy density for 
various time steps. Individual solitons as decay products are more apparent in this figure. 

The balance of H and V charges is of course respected in the decay (plj), which can 
be easily demonstrated by computing Qmv and Qmu for V^/, Hi, Vq , and H\. From the 
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FIG. 5: The decay of . (a) (\)\ and (\)2 as a function of x after the decay of . (b) The energy 
density as a function of x for various time steps. The dotted Une for Vq^ is added in order to clarify 
the trajectory of this decay product. 

conservation of energy point of view, M(yl) = 42.9, while M{Hf)+M{Vf^)+M{Hl) = 31.88 
which shows that the decay is allowed, and the excess energy is transferred to the kinetic 
energies of the solitons. As Fig. (i|(b)) shows, no appreciable energy is radiated away as 
small amplitude wavelets. 

Numerical results show that V2 is also unstable and decays via 

Vi Hi + Hi + + Hi + Hi (22) 

Figure (|^) illustrates the decay of ¥2- It can be seen from Fig. (|](b)) that a HIVqHI com- 
pound (i.e. Vi) is produced first. This compound subsequently decays into its components 
in a short time. Note that MlV-]-) = 89.6 and the sum of the rest energies of the decay 
products is Xlproduct ^ ~ 62.47. The difference goes to the kinetic energies of the decay 
products. The total energy and all topological charges are conserved. In order to demon- 
strate the conservation of topological charges, let us write down the charges of individual 
solitons: 

V2 : Qiv = +1 SI : Q2H = — 1 -^1 : QiH = —1 
: Q,v = +1 Hi : Q,H = +1 : Q2H = +1 
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FIG. 6: The decay of (a) cpi and 02 as a function of x after the decay of Vg*^. (b) The energy 
density as a function of x for various time steps. 

all other charges being zero. We thus have Qiv(LHS)=+l and (5iy(RHS)=+l, and all other 
charges equal to zero for both sides. 

Surprisingly enough, although the decay of via 



(23) 



is allowed by conservation of energy and charge, such a decay is not observed. We thus 
conclude that V-^ is a metastable soliton, and its fission needs some external trigger. 



V. CONCLUSION 



Although the Sine-Gordon equation is an integrable system, even minor modifications 
of the equation, usually exploit its integrability. We described in this paper and elsewhere 
0, that interesting and rich behavior may result by suitable modifications of the sine- 
Gordon equation irrespective of being integrable or not. Removal of mass degeneracy, soliton 
confinement, and soliton decay are among such properties. 

For the system introduced in this paper, we presented a class of exact, single soliton 
solutions in the particular case where the system reduced to the sine-Gordon equation (0i = 
or 02 = 0). Other classes of static solutions were computed numerically using a variational 
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algorithm. These static solutions were then fed into a numerical program which computed 
the dynamical evolution of the solution. We found that most of the static solutions were 
stable, while a few underwent decay into lower energy solitons. Topological charges where 
conserved throughout these dynamical processes, in addition to the conservation of energy 
and linear momentum which results from the invariance of Lagrangian under space and time 
translation. In the particular decay (|2^), two interesting effects were observed: First, the 
decay did not start promptly. Rather, V2 decayed into -ff^^^-^l first, and V-^ decayed via 
(PI) in a later stage. Second, the H2) pair apparently remained as a bound system in 
the decay of In order to check whether and H2 do form a bound system, we did 
numerical calculations (both variational and dynamical) for this system. It turned out that 
the system split into Hi and H2 which moved away from each other. We conclude that 
H1H2 do not form a bound system. 

Finally, let us discuss briefly the kinematics of solitons in the decay process. Since 
our system is relativistic, principles of conservation of energy and momentum should be 
applicable in their relativistic form. Ignoring the energy and momentum radiated away in 
the form of small amplitude wavelets, we may write 

n 

M^ = Y,1^^^C\ (24) 

i=l 

and 

n 

J^7imit;i = 0, (25) 

i=l 

for the decay 

S ^ Si + S2 + ... + Sn. (26) 

In equations (p^) and (p5|), M is the mass of the S soliton, are the masses of Si solitons, 
Vi are the corresponding velocities, and 7, = (1 —vf)^^^"^. The system of equations (^) and 



(p5|) have a unique solution for fj, if n = 2. For n > 2, we have only two equations for n 
unknowns and the equations do not have a unique solution. However, the decay pattern and 
the distribution of velocities among the decay products seem to be predetermined in our 
numerical results (compare the decay pattern of in Figures (^ and (^)). If we note that 
the decay to the daughter solitons does not happen in a single stage, but rather it proceeds 
in successive stages, it becomes clear why the distribution of velocities among the decay 
products follow a unique pattern. In each stage, an unstable soliton decays into two decay 
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products which is actually observed numerically. The velocities of the two decay products 
is unique, according to ( p^ and (pS]). One of the decay products, in turn, decays into two 
solitons in later stage, and so on. 
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